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1. Introduction 

The striking observations of large elliptic flow and jet quenching at RHIC [|I|] hints at an early 
equilibration of the dense matter formed in relativistic heavy ion collisions. This in turn suggests 
that the constituents of this medium, which at RHIC equilibrates at temperatures less than about 
300 MeV or energy densities of about 15 GeV/fm^, are still interacting strongly. The equation 
of state (EoS) describing this medium thus should exhibit significant deviations from ideal gas 
behavior. This has indeed been seen in calculations of the QCD equation of state on the lattice, 
where deviations from the non-interacting limiting behavior, e = 3p, show up as a large peak in the 
trace anomaly in units of T'^, &''^{T)/T'^ = {£ - 'ip)/T^- 

Pretty soon bulk thermodynamics, and in particular the EoS, will be probed at even higher 
temperatures at the LHC. This will be complemented by experimental studies of matter at lower 
temperature but non-zero baryon number density through low energy heavy ion collisions at RHIC 
and at the future European heavy ion facility, FAIR, at GSI/Darmstadt. It thus will be of importance 
to get control over the EoS in a wide range of temperatures and to analyze also its dependence on 
the baryon number density or the baryon chemical potential. 

While good control over the continuum extrapolation of numerical results on the EoS has been 
reached in lattice calculations for the SU(3) gauge theory, much work is still needed to get similarly 
good control over calculations performed with physical light and strange quark masses. We will 
discuss here the current status of these lattice calculations. 

During recent years the steady progress made in lattice calculations through important new 
algorithmic developments [g|] and increased computing resources helped to greatly advance our 
knowledge of the thermodynamics of strongly interacting elementary particles [||]. On the one 
hand, lattice calculations of the equation of state [Q ||, ^ and transition temperature [0, [s], ^ with 
an almost realistic quark mass spectrum are now possible at vanishing chemical potential. On the 
other hand, lattice calculations based on various approximation schemes now also provide first 
insight into the thermodynamics at non-zero quark or baryon chemical potential {}JLq) [[To|]. 

In this review we will discuss some of the results obtained recently with ^(a^) improved stag- 
gered fermion formulations. As these results have been presented in July /August of 2007 in a very 
similar format at the '4th International Workshop on Critical Point and Onset of Deconfinement' 
and at the 'XXV International Symposium on Lattice Field Theory' the write-up of these talks has 
been splitted into two parts. In this first part we will focus on recent results on the equation of state 
at vanishing and non-vanishing chemical potential and discuss results on strangeness and baryon 



number fluctuations which may become detectable in heavy ion collisions. Part II [11] focuses 
on a discussion of results on the nature of the QCD transition and reports on the current status of 
calculations of the transition temperature in QCD. 

2. The QCD equation of state 

While continuum extrapolated results for the EoS of the SU(3) gauge theory have been estab- 
lished more than 10 years ago ||3|], this still has to be achieved for QCD with a physical quark mass 
spectrum. Lattice calculations with light dynamical quarks, of course, introduce their own set of 
problems into the numerical study of QCD thermodynamics. The well-known problems with the 
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Figure 1: Pressure in the SU(3) gauge theory calculated for the standard Wilson action (1 x 1) and several 
improved actions analytically in the infinite temperature, ideal gas limit on lattices with temporal extent A^^ 
(left). The right hand figure shows numerical results obtained on lattices with temporal extent A^^ = 4. 

discretization of fermion actions that require to break the chiral symmetry of the continuum action 
and/or requires the introduction of additional fermion species leads to a cut-off dependent modifi- 
cation of the hadron spectrum that will influence thermodynamic properties at low temperatures, 
i.e. in the hadronic phase. This introduces additional cut-off dependencies into lattice studies of 
QCD thermodynamics that are not present in pure gauge theories and require a detailed analysis. 
Aside from this, however, studies of QCD thermodynamics suffer from similar discretization er- 
rors as pure gauge theories which strongly influence studies of bulk thermodynamics in the plasma 
phase of QCD and, in particular, show up as a strong cut-off dependence of the high temperature 
limit of bulk thermodynamic observables. These cut-off effects essentially arise from the need to 
replace derivatives in the QCD Lagrangian by finite differences. In the next subsection we briefly 
recall the importance of these effects for the analysis of thermodynamics in the high temperature 
phase of SU(3) gauge theories before starting to discuss the thermodynamics of QCD with two 
degenerate, light quark masses and a heavier strange quark mass. 



2.1 Cut-off effects in QCD thermodynamics 
2.1.1 SU(3) gauge theories revisited 

In the standard approach to the discretization of a SU(3) gauge theory the Euclidean action is 
represented by a sum over (1 x l)-Wilson loops (plaquettes). This introduces &{(?') discretization 
errors. It is well known how to eliminate these leading order discretization errors on the tree level by 
adding larger Wilson loops, eg. (1 x 2) or (2 x 2) loops. This shifts the leading order discretization 
errors to ^(a^). Other improvement schemes, like for instance the Iwasaki action, improve certain 
non-perturbative features of the lattice action, but do not aim at an elimination of the &{a^) errors. 
These discretization errors also show up in lattice studies of thermodynamic quantities on lattices 
with finite temporal extent (A^t); as the relevant scale is the lattice cut-off in units of the temperature, 
dr = l/Nt, ff{a") cut-off errors thus show up in finite temperature calculations as ff{l/N") errors. 
This can explicitly been analyzed in the infinite temperature ideal gas limit of QCD. In Fig.|I](left) 
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Figure 2: Cut-off dependence of the infinite temperature ideal gas pressure in staggered fermion formula- 
tions of QCD thermodynamics. Shown are analytic results for the pressure calculated in the thermodynamic 
limit {Na -^ °°) on lattices with temporal extent A^^ for the standard 1-link action as well as two i^{a^) 
improved actions. Results are shown for vanishing and non-vanishing values of the chemical potential 



we show the cut-off dependence of the pressure calculated analytically for various SU(3) actions 
as function of A^^ and normalized to the continuum (A^^ — > °°) Stefan-Boltzmann value. 

One lesson learned in studies of the EoS of a pure gauge theory is that it is extremely important 
to use at least G{a^^ improved gauge actions for the calculation of bulk thermodynamic observ- 
ables like eg. energy density and pressure, if one wants to reduce cut-off effects at high temperature 
already on lattices with small temporal extent. As can be seen in Fig. |l|(right) also at temperature 
close to the deconfinement transition temperature cut-off effects on lattices with fixed temporal 
extent (A^^ = 4) follow the pattern seen already in the infinite temperature limit. While the pressure 
calculated with actions that have 0'{a^) cut-off errors deviates strongly from the continuum extrap- 
olated result also in the vicinity of the transition temperature calculations with iff{a^) improved 
gauge actions yield results on coarse latices that are already close to the continuum extrapolated 
result. 



2.1.2 Thermodynamics with staggered fermions 

These generic features also hold true in calculations performed in QCD with light quarks. The 
standard discretization scheme of the fermionic part of the QCD action (1-link actions) introduces 
ff{a^) errors that can be removed in thermodynamic observables on the tree level by introducing 
suitably chosen additional 3-link terms (Naik, p4). The resulting cut-off effects in the infinite 
temperature, ideal gas limit are shown in Fig. ^. We note that the improvement schemes improve 
bulk thermodynamic observables at vanishing as well as non-vanishing chemical potential. These 
actions are usually supplemented with additional terms that reduce the flavor symmetry breaking 
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(smearing) and/or reduce cut-off effects beyond the tree level (tadpole improvement) ^ None of 
these improvements, however, alter the behavior in the infinite temperature limit. They are thus not 
expected to be essential for studies of bulk thermodynamics in the high temperature phase of QCD. 



The early calculations with unimproved staggered fermions [ ]12[ ] and calculations with unim 



proved Wilson fermions [ |13[ ] as well as recent calculations with the l-link, stout smeared action 
action [Q] indeed showed discretization errors similar to those observed in pure gauge theories; at 
high temperature energy density and pressure deviate strongly from the continuum ideal gas be- 
havior and reflect the large cut-off dependence of the ideal quark-gluon gas in these discretization 
schemes^ . We therefore will focus in the following on a discussion of results on bulk thermody- 
namics obtained with iff{a^) improved gauge and fermion actions. 

2.2 QCD Equation of State at vanishing chemical potential 
2.2.1 Bulk thermodynamics: Trace anomaly and entropy density 

Lattice calculations of the QCD equation of state with almost physical light quark masses and a 
physical value of the strange quark mass have, so far, been performed only on lattices with temporal 
extent Nr = 4 and 6. These calculations have been performed on lines of constant physics (LCP), 
i.e. with bare quark masses chosen such that a set of hadron masses stays constant as the continuum 
limit is approached. In order to study QCD thermodynamics over a wide temperature range it is 
necessary to determine the parameters that characterize the LCP at zero temperature over a wide 
range of lattice cut-offs. On the LCP one, furthermore, has to calculate an additional observable 
that can be used to set the temperature scale for the finite-T calculations. As this scale has to be 
determined at zero temperature, i.e. on large lattices, for quite a few different cut-off values (a) 
corresponding to the temperature values (T = l/N-cCi) one wants to analyze, it is advantageous to 
use a gluonic observable, e.g. the static quark potential, that is easy to calculate, is not influenced 
strongly by the chiral sector of QCD and has a weak cut-off dependence. In this way it becomes 
possible to minimize the influence of spurious zero temperature cut-off effects on thermodynamic 
studies and to control the genuine cut-off and quark mass dependence of numerical calculations at 
finite temperature. 

In Fig. ^ we show results for the trace anomaly, @^^{T)/T'^ = (^ ~ 3/?)/r^, and the entropy 
density, s/T^ = {e + p)/T^, calculated with &{a^) improved gauge and staggered fermion actions 
[§, ^]. We note that the entropy density as well as energy density and pressure are not independent 
observables but are obtained from @^^{T)/T^ using standard thermodynamic relations, i.e. after 
having evaluated the trace anomaly on the lattice, one obtains the pressure through integration over 
the trace anomaly. 



r4 T,^-Jr/^^^- ^^-'^ 



'For details on staggered fermion actions used in large scale studies of QCD thermodynamics we refer to the 
literature on, eg. the l-link, stout smeared action [H], the asqtad action [pi] and the p4fat3 action |H] (and references 
therein). 

^It sometimes is attempted to correct for these cut-off effects by normalizing numerical results obtained at some 
value of the temperature in simulations on lattices with given temporal extent N-c to the Stefan-Boltzmann value for 
pressure and energy density on the same size lattices at infinite temperature. This, however, does not solve the problem 
as cut-off effects are known to be temperature dependent. Moreover, it leads to thermodynamic inconsistencies when 
calculating, for instance, the trace anomaly, (e — 'ip)/T , which vanishes in the infinite temperature, ideal gas limit. 
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Figure 3: The trace anomaly, &^^{T) = e — 3p, in units of T^ (left) and the entropy density, ,s- = e + p, 
in units of T^ calculated with the p4fat3 action [p| on lattices with temporal extent Nz ^ 4- and 6. For 
0/i/J/7"4 Yve also show results on A^^ = 8 lattices (diamonds) obtained at high temperature. For N^ = 6 
results from calculations with the asqtad action are also shown [@]. In the right hand figure we also show 
the temperature scale TrQ (upper x-axis) which has been obtained from an analysis of static quark potentials 
at zero temperature [^]. The MeV-scale shown on the lower x-axis has been extracted from this using 
ro = 0.469 fm. 



Choosing the initial temperature Tq deep enough in the hadronic phase where the pressure becomes 
exponentially small such that the contribution of p{Tq)/T^ can be ignored safely, one obtains the 
energy and entropy densities through suitable combinations of @^^ {T) and p{T). The small verti- 
cal bar in the right hand part of Fig. fright) shows an estimate of the uncertainty that arises from 
setting the pressure to zero at To = 100 MeV. It is based on an estimate of the entropy density in a 
hadron resonance gas at this temperature ||^. 

In Fig. ^left) we show results for the trace anomaly obtained in calculations with the p4fat3 
action on lattices with temporal extent N^ = A, 6 and 8 on large spatial volumes, i.e. V^I^T = 
Na/N-i = 4 or larger. The insertion in this figure shows the high temperature part of (e — l>p)/T^. 
Obviously cut-off effects are an issue in this regime. For r>400 MeV, i.e. at temperatures larger 
than about twice the transition temperature, results on the N^ = 4 lattice drop more rapidly than 
on the Nt = 6 lattice. The calculations on the A^^ = 8 lattice, however, confirm the latter results 
and suggest that cut-off effects are under control in this high temperature region already on the 
Nx = 6 lattices. This is consistent with the analysis of cut-off effects for the p4-action in the infinite 



temperature, ideal gas limit []14[]. 

Also at r ~ 200 MeV, i.e. at temperatures just above the transition region, the trace anomaly 
shows some cut-off dependence. The peak in (e — l>p)/T^ is much more pronounced in the Ni = A 
calculations than in the Nt = 6 case. This can be traced back to the rapid change of non-perturbative 
j8 -functions in the crossover region from weak to strong coupling ^. As can be seen from the 
entropy density shown in Fig. fright) this leads to a noticeable cut-off dependence in s/T^ in 
this temperature range. It, however, has little effect on the high temperature behavior. This also 
holds true for energy density and pressure [Q]. It also is reassuring that calculations on performed 
Nt: = 6 lattices with the p4fat3 action on quite large lattices, V^'^T = 4, are in good agreement 
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Figure 4: The ratio p/e as function of the fourth root of the energy density obtained from calculations with 
the p4fat3 action on lattices with temporal extent Nx=A and 6. Also shown is the velocity of sound extracted 
from a fit to p/e [0| and using Eq. 2.2 The dashed curve at low energies shows the results for p/e calculated 
in a hadron resonance gas model (HRG). 



with simulations that have been performed with the asqtad action on a smaller physical volume, 
yi/3 J _ 2 [^ Results from the latter calculation are also shown in Fig. ^ The asqtad action 
has quite a different cut-off dependence at high temperature, it uses non-perturbatively improved 
(tadpole) couplings and also incorporates a more sophisticated smearing of 1-link terms in the 
staggered action to reduce flavor symmetry breaking effects. The good agreement between asqtad 
and p4fat3 simulations thus suggests that these features only play a minor role in the common 
temperature range explored in both calculations, 150MeV<r<400MeV. 

The results shown in Fig. ^ have been obtained in calculations with a physical strange quark 
mass and light quark masses that are about (2-2.5) times larger than in nature. This difference is of 
no significance at high temperature as the quark masses are small in units of the temperature^ . It 
may, however, play a role in the low temperature hadronic phase. From the experience gained in 
simulations with different light quark masses [^, ^] it is to expected that the region of sudden rise 
in the trace anomaly as well as the entropy density shifts to somewhat smaller temperatures in the 
case of physical quark mass values. Cut-off effects will lead to a similar effect. This deserves a 
further careful analysis (see also disccusion in part II [pl[]). 

2.2.2 Equation of state and velocity of sound 

For the description of the expansion of dense matter created in heavy ion collisions, in partic- 
ular its hydrodynamic modeling, the temperature dependence of bulk thermodynamic observables 
is not of direct interest. It is more relevant to get good control over the dependence of the pressure 
on the energy density, p{e), and deduce from this the velocity of sound. 



Ap_ 
d£ 



Ap/e p 
de e 



(2.2) 



The renormalization group invariant light quark mass for the calculations performed with the p4fat3 action has 
been estimated to be m'^^' = 8.0(4) MeV [R]. 
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Figure 5: Light ixi) and strange (jj) quark number susceptibilities at vanishing chemical potential: The left 
hand side shows XijTT'^ calculated with the p4fat3 action on lattices with temporal extent A^^ = 4 and 6. Ver- 
tical lines indicate here the location of transition temperatures deduced from peaks in chiral susceptibilities 
[pp on lattices of same size. The right hand side shows Xs/T^ calculated with p4fat3 and asqtad actions on 
lattices with temporal extent A^^ == 4, 6 [ |l^ and 8. All results are from calculations with a physical strange 
quark mass and almost physical, degenerate light quark masses, mi = 0. Im,. Data for A^^ = 8 are preliminary 
data of the hotQCD collaboration [|l7[. 



This is shown in Fig. ^. We note that the velocity of sound gets close to the ideal gas value at 
energy densities of about 15 GeV/fm^ or temperatures T ~ 300 MeV. However, for smaller energy 
densities, i.e. in the entire energy range of interest to the current RHIC experiments equation of 
state, p{s), as well as the velocity of sound drop rapidly. In the transition region, i.e. at e ~ 
1 GeV/fm^ one finds c^ ~ p/e ~ 0.09. 

We thus conclude that in the entire density regime, relevant for the expansion of dense matter 
created at RHIC, the QCD EoS shows large deviations from the conformal limit, c^ = p/e = 1/3. 

2.3 Non-vanishing chemical potential 

At small values of the quark chemical potential corrections to the equation of state can be 



evaluated systematically in terms of a Taylor expansion [ ]15[ ] in light (/X/ = jU„ = Hd) and strange 
(Hs) quark chemical potentials. The expansion coefficients are evaluated at vanishing values for the 
chemical potential and are related to thermal fluctuations of light quark and strangeness numbers 
in a thermal medium at /i/ ^ = 0. The leading order, quadratic corrections are given in terms of the 
diagonal quark number susceptibilities Xi and Xs as well as an off-diagonal susceptibility Xis, 



p{T,Hi,n,) _ piT,0,0) ixi f^i 

7^4 ^ 2T^\T 
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(2.3) 



(2.4) 



In the infinite temperature limit these susceptibilities will approach in the continuum limit the value 
for massless two and one flavor quark number susceptibilities of an ideal quark gas, respectively, 
i.e. \miT^ooXi/T^ = 2, Xmvj^^Xs/T^^ = 1 and \miT^ooXh/T^ = 0. It is obvious from Fig. that 
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Taylor expansion coefficients of the pressure, eg. the quark number susceptibilities, suffer from 
similar cut-off effects as the pressure evaluated at /i/ ^ = and, on the other hand, profit also in the 
same way from a systematic ^{a^) improvement of the action. 

Preliminary results for light and strange quark number susceptibilities, Xs.h obtained with the 
p4fat3 action on lattices with temporal extent Nx = A and 6 [ jl6| ] and the asqtad action [^ are 
shown in Fig. |5[ The quark number susceptibilities change rapidly in the transition region. Like 
bulk thermodynamic observables, they are sensitive to deconfinement and reflect the rapid change 
from heavy hadronic degrees of freedom to light partonic degrees of freedom. In Fig. ^(right) we 
show in addition preliminary results for the strange quark number susceptibility obtained by the 
hotQCD collaboration from calculations on a lattice of temporal extent N^ = ^ with the asqtad 



and p4fat3 actions [17]. These results indicate a shift of the transition region to somewhat smaller 
temperatures, but otherwise follow the same pattern seen already in the calculations on the coarser 
A^T = 4, 6 lattices. We will discuss the deconfining features as well as the cut-off dependence seen 
in calculations of quark number susceptibilities in more detail in part II of this review. Here we 
only point out that ;^/ ^ have a temperature dependence very similar to that shown in Fig. ^ for 
the entropy density, which also is shared by energy density and pressure. For chemical potentials 
IJ-i/T < 1 they give the dominant finite-jU contribution to these thermodynamic observables, which 
stays below 10% for }Xi/T<Q.l. This clearly is different for higher order derivatives of the partition 
function, that are sensitive to thermal fluctuations at non-zero baryon number density. 

Quark number susceptibilities play an important role in the analysis of the QCD phase diagram 
at non-zero baryon chemical potential. If there exists a second order phase transition point at non- 
zero baryon chemical potential, the baryon (or light quark) number susceptibility will diverge at 
this point. It thus is important to get control over the dependence of X/,i on /i/. The leading order, 
quadratic corrections to Xi and Xs are given in terms of fourth order expansion coefficients of the 
pressure |Jl8|]; for vanishing }Xs one finds, 

-YT- = -Jt- + 12C40 [yJ + ^(^' ) ' ~Y2- = -J2- + 2^22 [y) + ^(M/ ) , (2-5) 

with expansion coefficients c„„, = (VT^)-^ [{dlnZ{T,Hi,n,)/d{ni/T)''d{njT)'") /nl/ml]^^^^. 

The expansion coefficients, C40 and C22, are positive and have a maximum at the temperature 
of the transition at vanishing /X/ ,. For non-zero light quark chemical potential the susceptibilities 
thus increase and start developing a peak that becomes more pronounced with increasing /X/. This 
is obvious from Fig. ^ where we show baryon number and strangeness susceptibilities rather than 
light and strange quark number susceptibilities. The former are obtained in analogy to the latter 
after replacing in the QCD partition function the {u,d,s) chemical potentials by appropriate linear 
combinations of {B,Q,S) chemical potentials [|16|, ^, ^. The leading order Taylor expansion 
for the baryon number susceptibility evaluated here as function of the baryon chemical potential, 
Hb = lUd + Ml, = 3/X/, for vanishing strange quark chemical potential, jJ-s = IJ-d — Ms. and also for 
vanishing charge chemical potential, /Xg = /!„ — /i^, is then given by ^, 

Xb{^b) = ^ iXii^B) + Xs(Mb) + 2Z/.v(Mb)) • (2-6) 



Other choices for the chemical potentials may be more suitable to resemble situations encountered, for instance, 
in heavy ion collision. Charge and strangeness chemical potentials may be adjusted to reproduce conditions met in a 
particular event-by-event analysis of fluctuations nl^. 
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Figure 6: Fluctuations of baryon number (left) and strangeness (right) at non-zero baryon chemical potential 
calculated in leading order Taylor expansion on lattices with temporal extent A^t = 4 in QCD with a physical 
strange quark mass and almost physical light quark masses, m/ = 0.1ms. 

We note that strangeness fluctuations do react to an increase in net baryon number (/x^ > 0). This 
also holds true when the strange quark chemical potential is not chosen to be zero but is tuned such 
to insure that overall strangeness vanishes [16, 21]. The preliminary results shown in Fig. ^ so 
far have only been obtained on lattices with temporal extent A^^ = 4. Of course, one has to check 
the cut-off dependence also in this case through calculations on larger lattices [p^. Moreover, 
truncation effects in the Taylor expansion need to be analyzed by going beyond leading order in 
the expansion. In fact, higher order corrections are important! Only in the next order does the 
position of the peak become sensitive to the chemical potential and will get shifted towards lower 
temperatures with in creasing /X/. 

3. Conclusions 

Using &{a^) improved staggered fermion actions with an almost physical quark mass spec- 
trum much progress has been made in calculating the equation of state of QCD and extracting 
phenomenologically important quantities like the velocity of sound and fluctuations of hadronic 
charges. The current studies, performed mainly on lattices with temporal extent A^^ = 4 and 6 show 
some cut-off dependence. Calculations performed at high temperature for the trace anomaly on 
lattices with temporal extent N-^ = S suggest that these effects are small. A more detailed analysis, 
in particular at smaller temperatures in the hadronic phase, however, is still needed before firm con- 
clusions on the EoS in the continuum limit can be drawn. Preliminary results obtained for quark 
number susceptibilities show that the transition region shifts to smaller temperatures with increas- 
ing N^. This, of course, in accordance with the shift seen in systematic studies of the transition 
temperature performed with the p4fat3 [^ and asqtad [^]. actions. To quantify the effect for the 
equation of state will require a more detailed study on lattices with temporal extent A^^ = 8. 
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